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We study the phase transition in SU(2)-Higgs model on the lattice using the 3D dimensionally reduced formal- 
ism. The 3D formalism enables us to obtain highly accurate Monte Carlo results, which we extrapolate both to the 
infinite volume and to the continuum limit. Our formalism also provides for a well-determined and unique way 
to relate the results to the perturbation theory. We measure the critical temperature, latent heat and interface 
tension for Higgs masses up to 70 GeV. 



1. WHY 3D SIMULATIONS? 

Perturbative calculations have been extremely 
succesful in describing the physics of Electroweak 
interactions at zero temperature. However, at 
finite temperatures a purely perturbative analy- 
sis fails because of infrared problems: it is well 
known that the effective potential of the scalar 
field cannot be computed perturbatively for small 
<f>, in the symmetric phase. Thus, the calculation 
of the quantities characterizing the phase transi- 
tion - for example, the critical temperature T c , 
interface tension a, and latent heat L - requires 
the use of non-perturbative methods. 

A direct way to include the non-perturbative 
effects is to perform 4D finite-temperature lattice 
simulations of SU(2)-Higgs models. However, in 
the interesting parameter range the theory is still 
weakly coupled, and we can use perturbative di- 
mensional reduction (DR) to convert the 4D ac- 
tion into a 3D effective one. This step consists of 
integrating out all the massive modes (not con- 
stant in imaginary time) of the theory. In this 
talk we present results from 3D simulations with 
Higgs masses up to 70 GeV (for earlier results, see 
[1,2]; the results presented here will be described 
in detail in [3]). 

We maintain that, in practice, 3D simulations 
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are the method of choice for studying the EW 
phase transition [4,5]: 

(I) 3D model has one or two essential mass scales 
less than the original 4D model: in 4D, the lattice 
spacing a and the linear size of the lattice N x have 
to satisfy the limits T <C a -1 <C mjj{T)N x . In 
3D, the heavy T-scale does not exist, and we have 
to require only that mw(T) < a" 1 < mjj{T)N x . 

(II) 3D theory is superrenormaltzable — this 
gives an exact relation between the 3D lattice and 
continuum couplings in the limit a — ► 0, and we 
can relate any lattice observable to the physical 
one for given Higgs and W masses. 

(III) For a given a and N x , the number of lattice 
variables is much less in 3D than in 4D, making 
the simulations easier. 

(IV) We can consistently include the effects of 
fermtons and even typical extensions of the Stan- 
dard Model (for example, minimal SUSY exten- 
sions, the two-Higgs model) to the purely bosonic 
3D SU(2)-Higgs simulations [5]. 

The dimensionally reduced 3D SU(2)-Higgs La- 
grangian is formally similar to the 4D one: 

L=\f 2 + (Di^(D^) + m|> 2 + A 3 (</> 2 ) 2 (1) 

where <f> 2 = <f>^<f> and the 3D couplings g% and 
A3 have dimension GeV (Here we discuss only 
the case where Ao — the temporal component of 
the gauge field — is integrated over). We relate 
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the 3D couplings to the 4D ones at 2-loop level 
by Green's function matching [4,5]; using this 
method the nonlocal 2-loop terms which plaque 
straightforward DR [6] do not appear at all. The 
Lagrangian (1) is an approximation of the exact 
3D one; by systematically estimating the effects 
of the neglected terms we can conclude that for 
mif>60GeV the errors are less than 1%, depend- 
ing on the observable. 

The 3D lattice action can be written as 

S = Pg ^(l-^Tri^) 

- PnY^^^-^+i (2) 

Due to superrenormalizability (II), we have an 
exact relation between lattice and continuum 
parameters (f3 G , f3 H , f3 R ) <-> (#|a, \ 3 /g'i, m\j g%) 
when a — ► 0; for example, (3cj = 4/(g|a) directly 
connects the coupling constant /3q to the lattice 
spacing a. In 4D, the corresponding relation con- 
tains the RG constant ALatt, which has to be 
fixed by measurements. The 3D parameters are 
parametrized as (h = m#/80.6GeV) [3,5] 

gl = 0.44015 T 

\ 3 / g 2 = -0.00550 + 0.12622 h 2 , . 

m\lg\ = 0.39818+ 0.15545 h 2 ( ' 

- 0.00190 h 4 - 2.58088 m 2 H /T 2 . 

Note that neither m# nor T here are true phys- 
ical quantities. The values m# = 35, 60 and 
70 GeV used here correspond to physical pole 
masses m H (T=0) = 29.1, 54.4 and 64.3 GeV in 
4D SU(2)-Higgs theory (without fermions). For 
some other 4D theory (but the same 3D one, see 
(IV) above) the physical masses would be differ- 
ent [3]. Due to their transparency and universal 
nature we discuss only the 3D-values in the rest 
of the paper. 

2. SIMULATIONS AND RESULTS 

We fix the parameters according to eqs. (3) and 
use rriH = 35, 60 and 70 GeV. For each m#, we 
use (3q = 5, 8, 12 and 20, which correspond to 
different lattice spacings. For each /3q we have 
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Figure 1. The distribution of R 2 = <f>^<f> for some 
rriH = 60 GeV, fta = 8 volumes. 

several volumes, allowing us to extrapolate the 
measurements (A) to the thermodynamical limit 
V —> oo and (B) to the continuum limit a — ► 0. 
All in all, we have 59 different combinations of 

m H ,l3 G , V. 

For each lattice we search for the transition 
point by adjusting We have mainly concen- 
trated our effort to the m# = 60 GeV case. In 
fig. 1 we show the distribution of R 2 = <f>^ <f> for 
the largest fta = 8 volumes at the critical cou- 
pling Ph,c- The first order nature of the transition 
is obvious. As a rule, our 3D results qualitatively 
agree with the 4D results [7-9] and the recent 3D 
simulation [10]. However, the statistical errors 
in 3D are considerably smaller. For reviews, see 
[11,12]. _ 

The critical temperature 

We monitor the phase transition with order 
parameters R 2 and L = J2 Xt i | Tr V xU x<i V x , 
where V is the SU(2) direction of the Higgs field 
<f> = RV . The critical coupling (3h,c is located 
with several different methods: 

(1) maximum of C(L) = ((L - (L)) 2 ) 

(2) maximum of C(R 2 ) = ((R 2 - (R 2 )) 2 } 

(3) minimum of the Binder cumulant of L: 
B(L) = l-(Li)/(3(L 2 ) 2 ) 

(4) Equal weight value for the distribution p(R 2 ) 

(5) Equal height value for the distribution p(L) 
For each individual volume, the definitions (1)- 
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Figure 2. The V —> oo limit of (3h,c measure- 
ments. 

(5) yield different values for (3h,c, but when V —> 
oo all converge to the same limit, as is shown 
in fig. 2. For other values of /3q the situation is 
similar. 

For each /3q, we convert the V = oo value 
of Ph,c to transition temperature T c . These are 
in turn extrapolated to the continuum limit, as 
shown in fig. 3. For m# = 60 we have high preci- 
sion data for /3q = 5,8,12 and 20, and a good fit 
requires that we use a quadratic fit in 1//?g- F° r 
rriH = 35 and 70 GeV linear fits are acceptable. 
The final results are given in table 1; in all cases 
the transition is unambiguosly of first order. 

Numerically, the T c values from the simulations 
are quite close to the perturbative ones, but due 

Table 1 

The critical temperature T c , the interface tension 
a and the latent heat L for different Higgs masses. 
The value of a at m# = 35 GeV comes only from 
/3q = 8 simulations. 



m ff /GeV 35 


60 


70 


T c /GeV 92.64(7) 
T c P ert /GeV 93.3 


138.38(5) 
140.1 


154.52(10) 
157.0 


cr/T 3 [0.0917(25)] 
a pevt/ T 3 o.061 


0.0023(5) 
0.008 


0.005 


L/T? 0.256(8) 

£ P erty T 4 o.22 


0.0406(7) 
0.041 


0.0273(16) 
0.028 
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Figure 3. The continuum limit of the critical 
temperature T c for m# = 60 GeV. Only the 
quadratic fit has an acceptable x 2 /d.o.f. 

to the very high accuracy, they still differ at ~ 
10(7 level, signaling significant non-perturbative 
and higher order perturbative effects. 

The interface tension and the latent heat 

We measure the interface tension with the his- 
togram method: at the critical temperature the 
distribution of the order parameter develops a 
double-peak structure (fig. 1). The interface ten- 
sion can be extracted from the limit 



T 1 Pm 

— = lim — log — 

T V^co 2A 6 P n 



(4) 



where A is the area of the interface and _P max 
and -P m ; n are the distribution maximum and the 
minimum between the peaks. To use eq. (4) finite 
size corrections are needed; for details, see [3]. 
A crucial requirement is the "flat minimum" in 
the distribution between the peaks; this excludes 
all but the largest cylindrical volumes from the 
analysis. 

In fig. 4 we show the V —> oo extrapolation 
of a for rriH = 60 GeV. These values are then 
further extrapolated to fta —> oo; the final value 
is a = 0.0023(5) T c 3 . This is substantially smaller 
than the perturbative result 0.008 T c 3 , and signals 
the presence of non-perturbative effects for a. 

For rriH = 35 GeV we cite only fta = 8 result 
(table 1), since we do not have "flat" histograms 
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Figure 4. The interface tension extrapolated to 
V —> oo for m# = 60 GeV. 
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Figure 5. The W mass around T c for m# = 
60 GeV and /3 G = 12. 



for other /3q values; the continuum value is quite 
likely considerably smaller. For m# = 70 GeV we 
cannot reliably extract any non-trivial value. 

The latent heat L can be extracted from the 
discontinuity of R 2 at T c . For details, we again 
refer to [3]; in contrast to a, the continuum limit 
can be taken for all mu, and the results are re- 
markably close to the perturbative values, as can 
be seen from table 1. 

The Higgs and W masses 

In order to measure m#(T) and mw(T) we 
perform a separate series of simulations around 
T c for rriH = 60 GeV. We observe a good scal- 
ing between /3q = 8 and 12. Both m#(T) and 
mw(T) have a discontinuity at T c , and the masses 
are higher in the symmetric phase. In fig. 5 we 
show mw{T) in units of g% = g 2 T. The value 
of mw(T > T c ) contradicts the analytical limit 
mw / ^3 < 0.29 [13]. Similar behaviour has been 
observed in 4D [9] and 3D [10] simulations at 
smaller mu- 

REFERENCES 

1. K. Kajantie, K. Rummukainen and M. Sha- 
poshnikov, Nucl. Phys. B 407 (1993) 356 

2. K. Farakos, K. Kajantie, K. Rummukainen 
and M. Shaposhnikov, Nucl. Phys. B 407 
(1993) 356 



3. K. Farakos, K. Kajantie, M. Laine, K. Rum- 
mukainen and M. Shaposhnikov, in prepara- 
tion 

4. K. Farakos, K. Kajantie, K. Rummukainen 
and M. Shaposhnikov, Nucl. Phys. B 425 
(1994) 67; Nucl. Phys. B 442 (1995) 317 

5. K. Kajantie, M. Laine, K. Rummukainen 
and M. Shaposhnikov, CERN-TH/95-226, 
IUHET-312, hep-ph/9508379 

6. A. Jakovac, hep-ph/9502313; U. Kerres, G. 
Mack and G. Palma, DESY 94-226 

7. B. Bunk, E.-M. Ilgenfritz, J. Kripfganz and A. 
Schiller, Phys. Lett. B284 (1992) 371; Nucl. 
Phys. B403 (1993) 453 

8. F. Csikor, Z. Fodor, J. Hein, K. Jansen, 
A. Jaster and I. Montvay, Phys. Lett. B334 
(1994) 405 

9. Z. Fodor, J. Hein, K. Jansen, A. Jaster and I. 
Montvay, Nucl. Phys. B439 (1995) 147 

10. E.-M. Ilgenfritz, J. Kripfganz, H. Perlt and 
A. Schiller, HU-BERLIN-IEP-95/7 

11. K. Kajantie, LATTICE '94, Nucl. Phys. B 
(Proc. Suppl.) 42 (1995) 103 

12. K. Jansen, this volume 

13. W. Buchmiiller and O. Philipsen, Nucl. Phys. 
B 443 (1995) 47 



